function l2 = L2(lam,l,phij)

global sigmaC sigmaL phi

l2 = l*(lam.^(-sigmaC/sigmaL))'.*(phij/phi(1)).^(1/sigmaL);